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Abstract 

We derive a massive analog of the ALD (Abraham, Lorentz and Dirac) equation, i.e., the 
equation of motion of a relativistic charged particle with a radiation reaction term induced 
by emissions of massive fields. We show that the radiation reaction term has a non-analytic 
behavior as a function of the mass M of the radiation field and both expansions with respect 
of M and 1/M are generally invalid. Hence the massive ALD equation cannot be written as a 
local equation with derivative expansions. We especially investigate the radiation reaction in 
three specific motions, uniform acceleration, a circular motion and a scattering process. 
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1 Introduction 



A charged particle emits radiation when it is accelerated and loses the kinetic energy. Effectively 
the emission induces a friction term in the equation of motion of the charged particle. The 
modified equation with the radiation reaction term becomes 



and it is called the ALD equation pp. The first term is the external force while the second one 
corresponds to the back-reaction of the radiation of electromagnetic fields. The coefficient a 
of the radiation reaction force is a = e 2 /6ir when the massless photon is radiated. The ALD 
equation itself is well established but it has an infamous problem of the runaway solutions. The 
ALD equation contains the third derivative term and the point particle is accelerated to the 
speed of light by the radiation force. One can remove such pathological solutions by imposing a 
regular boundary condition at the infinite future. But then, the solutions must be accelerated 
before the external force is applied. It is against the causality and called the problem of the 
preacceleration. A pragmatic resolution is to treat the backreaction term as a perturbation 
j2], and replace the ALD equation by the so called Landau Lifshitz equation. The runaway 
solution and the problem of preacceleration can then be removed but the validity of the Landau- 
Lifshitz equation is restricted to a situation when the acceleration is sufficiently small. It is also 
important to investigate how the correction to the ALD equation is obtained from the quantum 
field theory calculation in which there is no such pathological problems (see, e.g. [H HI [5j [6]). 

An important aspect of the ALD equation can also be inferred from the study of the radiation 
reaction in other situations such as in other space-time dimensions j7], in curved backgrounds 
or the backreaction due to the gravitational radiation [H [101 Ej- in ai l °f these situations the 
radiation reaction term becomes nonlocal, and a local differential equation such as the ALD 
equation can be obtained only in special limits. Hence it will be important to see whether such 
a nonlocal effect also arises when the radiation field becomes massive, and to see how such a 
nonlocal effect is related to the problem of runaway solutions. 

In this paper, we consider a relativistic point particle coupled with a massive scalar field, and 
investigate properties of radiation reaction by emissions of the massive scalar fields. This system 
is considered as a toy model of the ALD equation with massive photons in plasma. We derive 
an analog of the ALD equation, but the radiation reaction term becomes generally nonlocal 
unlike the massless ALD equation, and the derivative expansion (1/M expansion) can not be 
applied. Furthermore, we show that the massless limit M — > is nonanalytic: the coefficient 
of the radiation reaction term contains a logarithmic term like M 2 log M or M 4 log M. 

The paper is organized as follows. In section 2, we derive the massive analog of the ALD 




(1.1) 
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equation. We first develop several tools for the calculation and compare the radiation reaction 
in massless and massive cases. The backreaction (radiation reaction) term in the massive case 
behaves non-analytically as a function of the mass M and contains a logarithmic term like 
M 2 log M. We also show that the derivative expansion (namely an expansion with respect to 
1/M) is generally invalid. Since the term itself is not singular as a function of M, it shows 
that the radiation reaction term is essentially nonlocal. In section 3, in order to see the explicit 
behavior of the backreaction, we evaluate the radiation reaction term numerically for specific 
motions, uniform acceleration and a circular motion. We show that, in the uniform acceleration, 
the radiation reaction term behaves like M 2 logM near M — 0. We also discuss a scattering 
process (a motion when the external force is applied during a finite time interval) to see that 
the nonanalytic behavior is not specific to the uniformly accelerated motion. In appendix \M, 
we discuss the non-analytic behavior of the backreaction based on the properties of the Bessel 
function which appears in the propagator of the massive scalar fields. 



2 Massive ALD Equation and Radiation Reaction 

2.1 Derivation of the Radiation Reaction Force 

The action of a relativistic charged particle interacting with a scalar field is given by 
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S = -m J dr y^¥z^ + J d A x j(x;z)<f>(x) + J d 4 x-(d^(x)d^<p(x) - M 2 2 (x)) (2.1) 

where the scalar current j(x; z) is defined as 

j{x; z) = e j dr ^¥F^\x - z(r)). (2.2) 

In the paper, we set the speed of light c = 1. This model is considered as a toy model of 
a charged point particle in massive photon fields in plasma. The equations of motion of the 
coupled system of the position of the particle z(t) and the scalar field 4>(x) are given by 

77103" = -e(&> -zPjrd v -iP)<j>(z), (2.3) 
(d^ + m 2 )0(x) = j(x; z). (2.4) 

The parameter r is chosen to satisfy the gauge condition z^z^ = 1. The above equation of 
motion is consistent with this gauge condition. In order to derive the massive ALD equation, 
we write the scalar field as a sum of an external part and a self-interaction part, 4>(x) = 
(f>ext{x) + <f>seif(x)- Then the equation of motion (12.31) for z{r) becomes 

m F = F? xt + F» elf (2.5) 
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where the external force F£ xt is given by the external field 4> ex t{x) as 

FL{z) = -e(0" - z»z v d v - z^ ext {z). (2.6) 

On the contrary, the self-force (radiation reaction force) F se if is generated by the charged 
current of the particle itself, and written in terms of the induced field (f) se if(x) as 

Kif(z) = - z»z v d v - zn<Psei f (z). (2.7) 

Here the induced field <fi se if(x) is solved in terms of the current j(x'; z) by using the retarded 
Green function 



<f>sei f (x) = j d*x'G R (x : x')j(x'-z) (2.8) 

where G R (x, x') satisfies the equation 

(d fl d^ + M 2 )G R (x,x r ) = 5 {A) (x-x'). (2.9) 

Let us first discuss some general properties of the radiation reaction force without using the 
explicit form of G R (x, x'). 

Due to the Lorentz symmetry, the Green function can be written as a function of the distance 
cr = (x — X 1 ) 2 , 

G R (x, x') = 6{x° - x'°)G(a). (2.10) 
Substituting this expression in <j)seif{x), one can rewrite the self-force as 

,t(z°) 

Kifiz) = -e 2 (^ - z^d v - F) / dr'G R ((x - z(t'))% =z(t) , (2.11) 



where r(x°) is defined by the proper time r which satisfies z°(t) = x°. In taking derivatives 
with respect to x, one first need to evaluate the integral at general x, and then take the limit 
of a; to a point z(t) on the trajectory of the particle. But noting that G(a) depends on and 
only through their distance a, one can express the right hand side of ( 12.111) in terms of the 



x 



'/' 



quantities on the trajectory 



d^J T X dr' G{a)\ x=z{T) = J T dr' 2(z^t) - z^r'))^-G(a). (2.12) 

The integrand is a function of a = (z(t) — z(t')) 2 and its derivatives while the integral is 
performed over the proper time r' of the trajectory. Since a and s = t — r' are related to each 



3 



other, we can either express the integral (12. 12j) in terms of a or s. Either expression has its 
own advantage, so we will explain both expressions in the following. 

We first express the integral f !2.12j) and the self-force term (12. lip in terms of s = r — t' . 
First z(r) — z{t') can be expanded in a power series of s as 



y"(s;r) 



Z[T) — Z(T 



E 

n=l 



TV. 



2 3 

sz»{t)- 8 - F(r) + ^^(r) + 
2 u 



(2.13) 



where ^"^(r) = d n z p (t) / dr n . We fix r and change the variable from r' to s = r — r'. By 



using the gauge condition z p z^ 
check the following relations 



1, z p z p = and 'z p "z „ + z^Zn = 0, it is straightforward to 



<t(s;t 

da(s; r) 
ds 



s 2 (l 



—z^z u H ), 

12 M ; 

r.2 



2y p y p = 2s(l - -z% + 



Using the relation d/da = (ds/da)d/ds, one obtains 

y" d 



P self 



ds(P p 



z p (r))G(a(s;r)), 



(2.14) 



(2.15) 



ypy p ds 

where P p = 5 P — z p (t)z u (t) is a projection operator which satisfies P p z u (r) = 0. The surface 
term does not contribute because P p (y u /y p y p ) — > P p z v = near s = 0. So we can perform an 
integration by parts and get the expression of the self-force term in terms of s 



F M - e 2 



ds 



V 



G(a). 



ds \y p y P , 

In order to evaluate it, let us first expand the integrand in powers of s, 



(2.16) 



r self 



ds 



ds 



)} + 



G(a) 



s 



-z t * + ~(z» + z>*z p z ) + 



G(a) 



(2.17) 



As we will see later, for a massless scalar field where G(a) ~ 5(a), higher order terms of s 
vanish except for the mass renormalization and the radiation reaction term. Thus one obtains 
the ALD equation, which is written as a local equation. On the other hand, for a massive 
scalar field, G(a) only damps as an inverse-power of a at s — > oo. Then integrals of s for 
sufficiently higher orders diverge. Since the integral H2 . 161) is shown to be finite after the mass 
renormalization, the divergence of the coefficients for higher orders indicates that the derivative 
expansion of the radiation reaction term is not valid. Hence the massive analog of the ALD 
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equation cannot be written as a sum of local terms. We therefore need to evaluate the integral 
(12.161) directly without the derivative expansion, and for this purpose, it is more convenient to 
express the integral in terms of a instead of s. 

The Green function G(a) is generally a complicated function of a (e.g. see (I2.26P for a 
massive scalar field. Therefore it is sometimes more appropriate to express the integral in 
terms of a instead of s, especially in evaluating the radiation reaction term numerically. The 
integral (I2.16P is rewritten by changing the integration variable from s to a as 

ds 
da 

In the region < s < oo, the function a(s) is single valued and da/ds > is satisfied. Hence 
we can solve s as a function of a. By using the variable I = -Ja, the self-force can be expressed 
as 

ds 



P self 



da 



( ~2f^ 
da 



v \da 



G{a) 



(2.18) 



r self 



dl 



dl 



dl I dl 



dl 



G(l 2 ) 
G(l 2 ] 



(2.19) 



In the second equality, we have used the following relations 

ds > 2 
dl 



— = 1 + 



z v z v V 



o(l 3 



(2.20) 



The usefulness of the expression (12.191) becomes apparent when we consider, e.g. the renor- 
malization of the mass term m'z. The radiation reaction term, (12.161) or (12.181) . contains a 
divergence which is proportional to z^, and should be absorbed by the mass renormalization. 
For a massless scalar field with G(a) ~ 5(a) = S(s 2 ), it makes no difference whether we sub- 
tract the divergent term in ( 12. 16[) or (I2.18p . However, for a massive case, G(a) is a nontrivial 
function of a(z, z') = (z^(t) — ^(r')) 2 , and since a depends on z^ and t' in a complicated way, 
an integral of G(a) such as J dsf(s)'z^(T)G(a) depends not only on z^(r) but also on the details 
of the trajectory z(t'). On the other hand, an integral of G(a) such as J da f (a) z^ (r)G(a) is a 
constant multiplied by z^(t). Therefore the mass renormalization of the massive ALD equation 
should be evaluated not in the expression (12.161) but in (12. 18D . 



2.2 Massless Radiation 



Before considering a massive case, we briefly review the evaluation of the radiation reaction by 
a massless scalar field. For a massless scalar field, the retarded Green function is simply given 
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by 



G R (x-x') = 9(x° ~x' ) 5 -^. (2.21) 

2n 



By substituting it in (12.171) or (12. 19[) §, we can obtain the ALD equation 



e 2 



mz» = F» xt + — (z" + #\g%). (2.22) 



Here the divergent term 



oo 



I z^ (2.23) 



is absorbed by the mass renormalization. 

We would like to mention the derivation of the ALD equation from the quantum field 
theory calculations j5]. The authors evaluated two types of diagrams which contribute to the 
backreaction. One type of diagrams corresponds to the radiation processes. The contribution 
of this type to the backreaction turns out to be written in the form of ( 12. lip but with Gr 
replaced by G_ 

G-(x,x') = G R (x,x') - G A (x,x') = 6 ( x ° - x '°)- e ( x '° - x °) G ( a ^ (2 24) 

2 

where Ga{x,x') is the advanced Green function. With this replacement, the terms with even 
powers of s in (12. 17ft vanish and the remaining terms become finite. Another type of diagrams 
corresponds to the self-energy (including the mass renormalization) of the point particle. It 
has a form of ( 12. lip but with Gr replaced by G + 

G + (x, x') = G R (x, x') + G A (x, x'). (2.25) 

It gives a divergent term which can be absorbed by the mass counter term. Summing these 
two types of diagrams, one can reproduce the ALD equation §. 



*The integration seems subtle since the delta function is just at the edge of the integration range. But the 
subtlety can be removed by going back to the original equation (|2.1ip . Here the integral is performed over the 
region t(x ) > r' > — oo and the delta function is always inside the region. 

§One might expect to subtract a divergence by replacing Gr by G_ or Gf also for the massive case discussed 
in the next subsection. But the Green function in the massless case depends on a only through 5(a), and the 
acausal time dependence of G a inside G_ and Gf does not appear in the final results. However, in a massive 
case, the Green function has a tail and a simple replacing of Gr by another type of Green functions such as G± 
violates the causality explicitly. 
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2.3 Massive Radiation 



The retarded Green function of a massive scalar field is given by 




9(x° - x'°) 



25(a) - 9(a) 



MJ x (M^a~) 
7^ 



(2.26) 



The first term in the parentheses is the same as the massless case while the second term 
represents a modification of the Green function by the mass of the radiation field. Since the 
second term disappears in the massless limit, one might expect that the massive ALD equation 
can be written as a perturbation of the massless ALD equation by higher derivative terms. 
The situation is not so simple, however. By the dimensional analysis, the coefficient of the 
derivative expansion must contain an inverse power of M, but such an expansion is shown to 
be invalid because of the nonanalyticity at M = 0. 

To see such nonanalytic behaviors of the derivative expansion, let us try to evaluate each 
coefficient of the derivative expansion by substituting the Green function (12.261) in (I2.19p . The 
leading order term is proportional to which corresponds to the mass renormalization. The 
next order term gives the coefficient of the ALD term. The coefficient, however, vanishes due 
to the cancellation between the massless propagator 5(a) and the massive modification of the 
Green function, 



This cancellation is universal and does not depend on the value of M. We are tempted to 
conclude that the ALD term vanishes for a massive scalar field, but it is not correct since the 
derivative expansion is shown to be invalid. It suggests that something strange happens in the 
massless limit. Let us evaluate the coefficients of higher order terms. They are written as 



Since the Bessel function behaves as J\(x) ~ l/y/x for large x (see Appendix A), the integral 
diverges for n > 1/2. Therefore the derivative expansion (or equivalently a power-series expan- 
sion with respect to I = ^/a) is not justified. One can show that the expansion in terms of s in 
(12.171) is also divergent. Such singular behavior of the derivative expansion in the massive case 
shows a peculiarity of the massless ALD equation where the radiation reaction force is written 
local derivative term. 

Though one can not write the radiation reaction term in terms of local derivative terms, 
the backreaction itself is finite even in the massive case. To confirm the finiteness of the 
backreaction, we go back to (12.181) or (12.161) . The integration can be divided into two parts, one 
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(2.27) 





involving the integration of 6(a) and the other involving the integration of Ji(M^fa). Since the 
first term is the same as the massless case, it is sufficient to show the finiteness of the integrals 
of the second part, 



oo 







and 



da ^MM^) 
da \ a 



Note that for general trajectories P^(d/da)(y u ds/da) or ds/da is regular at finite a and falls 
off faster than a -1 / 2 at a — > oo(or s — > oo). On the other hand, M Ji(M \fa) / \fa is regular at 
finite a and behaves like 

-cos^ (2.31) 
a 

at infinity. So the integration (12.181) is finite for a fixed M, and we get a finite amount of 
backreaction by emissions of a massive scalar field as expected. We see this finiteness explicitly 
in the next section for special types of trajectories. 

We have seen that, if we perform the derivative expansion, the coefficient of the ALD term 
vanishes for a finite M while it becomes (e 2 /127r) at M = 0. One may also be interested in the 
mass expansion near M = of the backreaction term and how the local equation can be derived 
from the massive ALD equation with a nonlocal radiation reaction term, (12. 18[) or (I2.16p . But 
it is not so simple as one might expect. For example, consider the convergent integral (12.291) . 
If we naively expand it as a power series of M, we get 

M L p ^(^)g aV+2) ■ (2 ' 32) 

But the higher order terms of M are divergent, and such an expansion is not valid. The 
backreaction must have a non- analytic behavior at M — 0. We will investigate this kind of 
non-analyticity in the next section. 

2.4 Mass Renormalization 

In this subsection, we discuss the mass renormalization. Because of the non-locality of the 
backreaction term, the mass renormalization becomes more subtle than the massless case. Since 
the non-local term is convergent, one may naively renormalize the mass only by the divergent 
integral. But then the renormalization of mass becomes independent of M and survives even 
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in the infinite M limit where the emission of such a massive field is highly suppressed. More 
explicitly, because of the identity of the Bessel function, 

lim M6{a) Jl ^ M ^ = 25(a), (2.33) 



the Green function vanishes in the large M limit 

llm Gr(X: x <) = llm 'fr" -*") U*) - e(a) MJ ^) = 0. (2.34) 

The backreaction therefore must vanish in the limit M — > oo. Thus we also need to take 
into account the non-local but convergent integrals, ( 12.291) and (12.301) . in discussing the mass 
renormalization at finite M. 

In order to make a consistent mass renormalization, we first introduce a function /(/) and 
define the renormalized mass term by 

POO 

,2 / ji ■■Uf(i\nn2 



5m z" = -e 2 / dl z"f(l)G(l 2 ). (2.35) 
Jo 

Then the self-force term can be divided into the mass renormalization 5mz^ and the radiation 
reaction force Ffi 0jSe if as 

Keif = Ffaj + SmzP. (2.36) 

The function f(l) must satisfy the condition f(0) = 1/2 so that the divergent term is correctly 
subtracted in (12.161) . It is consistent with the mass renormalization in the massless case. For 
the massive scalar field, the above subtraction (12.351) has a contribution not only from the 
term 5(a) in the Green function but also from the non-local term ~ J\ (My/a) j \fa . Another 
condition for /(/) is that the mass renormalization should vanish in the large M limit. The 
simplest choice is /(/) = 1/2, but the choice is not unique. The ambiguous part of the mass 
renormalization is written as 

Sm ^^r m) Mim. ( , 37) 

Jo ' 

It is independent of the trajectory of the particle. Hence a different choice of /(/) does not affect 
how the backreaction term depends on the details of the trajectory such as the acceleration or 
the velocity of the particle's motion. In this paper, we take the simplest choice f(l) = 1/2. 
In the next section, we will show, if we renormalize the mass by /(/) = 1/2, the radiation 
reaction term becomes o(a 3 ) when a particle is uniformly accelerated with an acceleration a. 
For other choices, it is generally proportional to o(a). This seems to justify the simplest choice 
f(l) = 1/2- 
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3 Radiation Reaction for Specific Trajectories 



In order to investigate the non-analytic properties of the backreaction, we evaluate ( I2.16P and 
(12.181) explicitly for some specific trajectories. We consider three uniform acceleration, 

a circular motion and a scattering process. 

3.1 Uniform Acceleration 

A uniformly accelerated point particle is described by the trajectory 

= (- sinh ar, - cosh ar, 0, 0), (3.1) 
a a 

where a is the acceleration and z v z v = —a 2 . i M (r) and z^(r) form a two dimensional vector 
space and other derivatives are written in terms of them as follows 

z (2n)„ = a 2(n-l)«M z (2»+l)/i = = q (3^ 

For the trajectory (13 .11) . the radiation reaction term vanishes 

Z" + ^ = 0. (3.3) 
l(s) or a(s) has the following simple form 

2 ^ a(r - r') 



Z = v^= ^(^(r) -^(t'))(^(t) - ^(r')) = -sinh v - (3.4) 
The backreaction becomes proportional to z^ since P^z v = and is given by 

£M f 00 / l + a 2 Z 2 /2 A 
F£„ = —F(a,M) = e 2 a(sinh ar, cosh ar, 0,0) / dZ 1 ] G R (l). (3.5) 

Because of (13.31) the delta function term in ( 12. 26ft vanishes in Ff^f. Then F(a, M) becomes 
r,, „^ 2 P ( l + a 2 l 2 /2 l\ MJAMl) 



2(^1 + gt 2 ) 3 2 J 4vrf 
= e 2 aMF(g), (3.6) 

where g = a 2 /AM 2 . In the second equality, we changed the integration variable from Z to 
t = Ml. All information of the backreaction to the uniformly accelerated particle is contained 
in the function F(a, M) or equivalently in F{g). 
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Figure 1: The function F(g) is plotted as a function of g = a 2 /AM 2 . 



The function F(g) can be expressed by using the Meijer's G-function (see for example |1 1 j ^ 

1 



F(g) 



i 



57T 



o {^/l + gT 2 ? t 



1 

8tt Stt^ 



g r 21 / J_|l 
7T 3 / 2 13 I 4^1 3/2,3/2,1/2 J * 



(3.7) 



The behavior of the function -F(g) is plotted in figure. [TJ In the two limiting cases of g = 
and — > oo, it can be approximated as 



limFfo) 



1 



167T 



<7 + 



(3.8) 



Substituting g = a 2 /AM 2 in F(g), one obtains the dependence of the backreaction on the 
acceleration a of the point particle and the mass M of the radiation field. 

Let us first look at the behavior of the backreaction F(a, M) as a function of M with the 
acceleration a fixed. It is plotted in Figure. El From ( 13.81) . one can see that near M = the 
backreaction F(a, M) becomes 

- 2 , e 2 aM e 2 M 2 



lim F(a, M) = lim e 2 aMF(— — 



S7T 



47T 



log M + ■ ■ 



(3.9) 



Note that e 2 aMF(g) vanishes at M = 0. However it is not analytic at M = and contains a 
logarithmic term proportional to M 2 log M. This is the reason why one could not expand the 
backreaction term with respect to M. One can also obtain the large mass limit, M — > oo, as 

1 e 2 a 3 



lim F(a,M) 

M— >oo 



16tt AM 



+ ■ 



(3.10) 
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Figure 2: aMF(g) as a function of M, with a = 1 



2.0 2.5 



Figure 3: aMF(g) as a function of a, with M = 1 

The backreaction vanishes in the limit as expected. 

We then fix M and look at the dependence of the backreaction on the acceleration a. The 
behavior of F(a, M) = aMF(g) as a function of a is plotted in Figure. [31 In the limit a — > oo, 
the back reaction term -F(a, M) becomes 



e 2 aM 

lim F(a, M) = lim ( h e 2 M 2 log a). 



Stt 



(3.11) 



It is proportional to a in the leading order, and has a correction of log a. The logarithmic factor 
log a has the same origin as the logarithmic factor logM in (13.91) . In the limit a — > which 
corresponds to the limit of the point charge staying at the origin, the backreaction is suppressed 



as a" 



e 2 a 3 



lim e 2 aMF(q) = — lim 



->> 



(3.12) 



If we took a different choice of f(l), the coefficient of the term proportional to a did not vanish. 
This implies the plausibility of the simplest choice /(/) = 1/2. Indeed, if we choose a different 



' It is written by an inverse power of M, but it does not mean that a derivative expansion is valid since an 
explicit trajectory was used to derive the result. 
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f(l), the change of the backreaction term is given by 

e 2 aM J dt (f(t/M) - = e 2 aSF(M). (3.13) 

It is proportional to a multiplied by a function of M. Though a different choice of f(l) changes 
the mass renormalization and the leading behavior of the radiation reaction in the a — > or 
a — > oo limit, it does not change the logarithmic behavior (log a or logM). 

Another important feature is the sign of the backreaction F(a,M). As shown in figure 
[21 F(a = 1,M) becomes either negative or positive, and such a behavior of F(a,M) is not 
special to f(l) = 1/2, but also holds for a general choice of /(/) (see the appendix iBl) . Since 
the backreaction is given by Ff^f = (z M /a)F(a, M), a negative value of F(a,M) is physically 
natural since it means that the backreaction term suppresses a change of the particle motion. 

On the contrary, if F(a, M) is positive, it enhances the acceleration when the particle is 
accelerated. Hence it is an indication of an instability of particle motion or a breakdown of 
the framework which we are using. The positiveness F(a, M) > for a specific region of the 
parameters is generally satisfied for a 3> M. It is interesting to investigate a relation with the 
problem of runaway solutions in the massless ALD equation. 

3.2 A Scattering Process 

In the previous section, we have seen a non-analytic behavior of the backreaction for a uniform 
acceleration. In this case, the point charge is accelerated for an infinitely long time, and the 
total radiation becomes infinite. Hence one may suspect that such a non-analytic behavior is 
caused by the infinite acceleration. In this section, we see that a similar non-analytic behavior 
does arise even when a particle is accelerated for a finite time interval. 

Now we consider a charged particle with a constant four velocity t> M at r < 0. During a 
finite time interval < r < L, the particle is scattered (or feels an external force) and then 
after r > L it goes away with a velocity v ' M . The trajectory is given by 

v M r, for t < 



for < t < L (3-14) 



q-[T) 

v'^r + c M , for r > L 

where q^{r) is a function satisfying the boundary conditions ^(0) = and q^(L) = c M + v ,fl L. 

At t < 0, the backreaction is since no radiation has been emitted yet. It is confirmed 
from the equation (12.161) . Since = sv^, the backreaction trivially vanishes by the projection 
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operator Pj^ = 5^ — z^iy. During the time interval < r < L, the point particle is scattered by 
the target (or source of the external field). Let us divide the integral (12.161) into two parts, one 
with ds and the other with ds. The first part does not cause any non- analytic behavior 
since the integration range J Q r ds is finite 0. So we are interested in the second part, 

? f°° , ( „ d y u , A MJAMl) 

F ™ = ~ e L ds v'^h + z{T> ) -^tr 1 - (3 - 15) 

Using the following relations 

= z*(t) - z"(t') = q"(r) - uV, = = v", 

ds 

° = y"y, = ?% + r' 2 - 2g%r', (3.16) 
one can rewrite the integral as 

e2 r ds m - - \ mj i{ mi) 

Using the property of I ~ r' ~ s at s — >■ oo, the integration becomes 

F£^-e 2 J°°dl (V + (ff - ^)Kg a - gV> a i(l - 2 ^ + — )) ■ (3.18) 

The leading power of / at / — >• oo in the parenthesis is given by the term proportional to l~ 3 . 
Due to the property of an integral of the Bessel function (see Appendix 13} , the integration 
gives a term proportional to M 4 logM. This term vanishes in the massless limit, but it is 
non-analytic and we cannot expand the backreaction term with respect to the mass M. Hence, 
even for a scattering process where the point particle is accelerated only during a finite time 
interval, the backreaction behaves non- analytically at M — 0. 

3.3 A Circular Motion 

The third example we consider is a circular motion. When the particle is moving on a circle, 
the trajectory is given by 

z^ = (77-, pcos^ujT, psin-fUT, 0). (3.19) 

The parameter p is the radius of the circular motion and pu(< 1) is the velocity of the point 
charge. These parameters (p, u, 7) must satisfy 

7 = 7==^ ' (3.20) 
y 1 — p z u 2 



" Since the integral is always convergent, we can obtain a power series of M n by expanding MJ\{Ml) with 
respect to M. 
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in order to hold the gauge condition = 1. Since i M , z M and z M are all independent of each 
other, the backreaction can be generally written in the form 

Fself( Z ) — e2 [-^n(0, COS 'JUT, Sm'JUT, 0) + F ALD (pU, - SmjUT, COS JUT, 0)], (3.21) 

where the first term is proportional to — z M and the second term is proportional to — (z^ + 
z^'z v z v ). A term proportional to i M vanishes due to the projection operator Pjf. The following 
relations are useful 

"z^ + z^'z v z v = pj 5 u 3 (— pu, sin7a;r, — COS7WT, 0), z v z v = —p 2 j 4 u 4 , 

a = {z(tY - z(t'Y)(z(t), - z(r%) = 7 V - 2p 2 (l - cos 7 ^). (3.22) 

The coefficient F m is given by 

f°° , f (1 + p 2 u 2 ) (cos 7 us- 1)+ 7 us sin 7 us 2 2 ,1 dl A MJAMl) 

F m = - dsp[ — : v +7W (-— - 1) — , 

J \ (7s — p z u siwyusj^ 2 as J Attl 

(3.23) 

and determines the mass renormalization. The other term Fald gives a coefficient of the ALD 
term 

f°° jus cos jus — sin jus MJ\ (Ml) pj^u 3 

Fald = - ds p— — : — — — . 3.24 

J (js — p z u sm jus ) z Airl llix 

The second term — (pj 5 u 3 /12ir) is equal to the ALD force when the radiation field is massless 
while the first term gives a correction in the massive case and comes from the non-local part 
of the Green function. 

In order to see the dependence of the backreaction on M, we first note that / ~ 7s at 
s —> 00. Then both of the integrands of Fald and F m are proportional to sinul or cos ul at 
/ — > 00 lj. Thus there are no non-analytic terms such as log M(see appendix |X] for details) in 
the case of the circular motion. 

Finally let us see how the radiation reaction changes as a function of M. In Figure. HI the 
ratio of Fald to the radiation reaction in the massless case is plotted. It becomes 1 at M — > 
and decreases as the mass of the radiation field M increases. Eventually it vanishes at M — > 00. 
This result is consistent with our expectation that the radiation reaction is suppressed by the 
effect of the mass of the radiation field. In Figure. [5] we plot the same ratio as a function of 
u. The ratio becomes 1 in the relativistic limit pu — > 1. When the frequency u decreases and 
the particle's motion becomes nonrelativistic, the radiation reaction is more suppressed in the 
massive case. This can be naturally expected since in the nonrelativistic region with a small u 
the emission is highly suppressed by the mass of the radiation field. 



** There is also a non-oscillating term ~ I 2 in F m . But as we see in the appendix [3] such a term with an 
even power of / does not produce a logarithmic behavior at M = 0. 
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Figure 4: F ALD /F^ sless = F ald /(-P7 5 w 3 /12tt) as a function of M, with p = 1 and u = 0.8 

0.8 - / 



0.6 - 



Figure 5: F ALD j 'F™g^ less as a function of tu with p = 1, M = 2 

4 Conclusions 

In this paper, we investigated the radiation reaction of a charged particle interacting with 
a massive scalar field. We first obtained a massive analog of the ALD equation. The most 
important observation is that the massive ALD equation cannot be written as a local equation 
with higher derivative terms. The coefficients of higher derivative terms are divergent and the 
derivative expansion of the radiation reaction term is invalid. Only in the massless limit the 
ALD equation becomes local. This sounds strange since the massless limit is more sensitive to 
the infrared effect and the nonlocal effect seem to become more important than a massive case. 
Technically speaking the locality arises in the massless ALD equation because the retarded 
Green function of a massless field has a support only on the light cone and is proportional to 
£(a"). On the contrary, the Green function of a massive field is distributed around the light 
cone and thus the radiation reaction term becomes nonlocal which is written only as an integral 
form. 

We also studied the nonanalytic behavior of the radiation reaction term in various situations. 
First we showed that it has a non-analytic behavior as M p logM 2 at M — > 0, where p is a 
positive integer depending on the details of the trajectory of the particle. Such non-analytic 
behavior generally appears at M — > even when the particle is accelerated during a finite 
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time interval. In appendix |A] we studied how such a non-analytic behavior with logM appears 
by looking at various integrals of the Bessel function. We also evaluated the backreaction in 
specific motions, a uniform acceleration and a circular motion. In both cases, we evaluated the 
backreaction (radiation reaction) as a function of the mass M and showed that it is suppressed 
when the mass becomes large. It is consistent with our physical intuition that the emission is 
suppressed in the massive case and accordingly the backreaction is reduced. 

Finally we would like to comment on a possible resolution to the runaway solutions in 
the massless ALD equation. As we saw at the end of Section 13.14 the positive region of the 
backreaction F(a, M) suggests that an instability will occur in the region a 3> M. It may 
be related to the pathological behavior of the massless ALD equation. In the massive case, 
when a particle is accelerated from a nonrelativistic region, the backreaction term F(a, M) is 
negative and suppresses the acceleration. Hence it tends to avoid the runaway type solution. 
It is interesting to investigate that if one can avoid the pathology of the runaway solution by 
introducing a small mass of the radiation field and then taking a massless limit. 

Acknowledgments 

We would like to thank S.V. Bulanov and K. Seto for discussions. The research by S.I. is 
supported in part by Grant-in- Aid for Scientific Research (19540316) from MEXT, Japan. We 
are also supported in part by "the Center for the Promotion of Integrated Sciences (CPIS) " 
of Sokendai. 

A Some properties of J\{x) 

The Bessel function J\(x) is given by the series 



which is an odd function of x. It is regular on the whole complex plane except at infinity. The 
asymptotic behavior of J\{x) is given by 



Though it falls off very slowly as at infinity, the oscillating factor makes the following 

integral convergent 






(A.2) 




(A.3) 
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In this appendix, we focus on the behavior of the function 



m i) = MJ fl (A.4) 
and its integrals in the two limiting situations, M — > oo and M — > 0. 

A.l M — > oo limit 

We first show that the function I(M, I) has the following property 

lim I(M,l) = 25{l 2 ). (A.5) 

M— >oo 

In order to show this, we introduce a continuous and regular function g(l) in the region < 
I < oo with the boundary condition linn^oo g(l) = 0. Then the following integral becomes 



poo POO 

lim / d(l 2 ) I(M,l)g(l 2 ) = lim / dl 2MJ 1 (Ml)g(l 2 ) 

I^oo J Q M-+oo J Q 

POO 

= lim 2 / dsJ 1 (s)g(s 2 /M 2 ) 

A/^oo J Q 

POO 

= 2 / dsJtisMO) = 2(7(0). (A.6) 



o 



In the last equality we changed the ordering of limM^oo and / °° ds. If it is justified, the equality 
(1A.5|) is proved. To justify it, one needs to show that 



lim J 1 {s)g{s 2 /M 2 ) = Ji(a)^(0) 

M-»oo 

is uniformly convergent at < s < oo. Namely we need to show that for any e there exists N t 
such that the inequality 

|Ji( S ){^ 2 /iV e 2 )-^(0)}|<e, (A.7) 

is satisfied for all < s < oo. 

Such N e can be found as follows. Let us denote the maximum value of \g(s) — g(0)\ for 
< s < oo by 5g max . Then, since Ji(s) behaves like 1/y/s at large s, there exists S e such that 
\Ji(s)Sg max \ < e is satisfied for s > S e . On the other hand, for the region s < S e , we can always 
find N e such that \g(s 2 /Nf^J — g(0)\ < e since g(l) is continuous at / = 0. To summarize, for 
the region s > S e , \Ji(s){g(s 2 /N 2 ) — g(0)}\ < \ Ji(s)Sg max \ < e is satisfied and for the region 



s < S e , | Ji(s){(g(s 2 /N 2 ) -g(0)}\ < \g(s 2 /N 2 )-g(0)\ < e. Thus the uniform convergence (jA7f]l 
is proved. 
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A. 2 Massless limit: M -» 

We then investigate the behavior of the integration 

/■oo 

Q(M)= / dlI(MJ)g(l), (A.8) 

in the massless limit M — > 0. One can similarly prove that 

hm /(M, /)<?(/) = (A.9) 

is uniformly convergent at < I < oo. We can exchange the ordering of the limit and the 
integral and show that 

Q(0)= hmQ(M) = 0. (A.10) 

Hence the non-local term like (I2.29P vanishes in the massless limit and it is reduced to the 
massless ALD equation. 

However, the analytic behavior of Q{M) is not so simple since logarithmic terms like 
M 2 logM are expected to appear around M — 0. In the following, we show how to derive 
such a non-analytic term of Q(M). First note that the integrand I(M,l) is analytic for both 
parameters, I and M. Then an integration over a finite range 

Q{M-L)= [ dl I(M,l)g{l), (All) 
Jo 

is also an analytic function of M . Non-analytic behavior like log M appears only when we take 
the integral region infinite A — > oo. 

Let g{T) be an analytic function at finite / which falls off as l~ n at infinity. Here n is a 
non-negative integer. We encountered two examples in the analysis of the radiation reaction. 
In the case of the uniform acceleration (13.61) . an integral with g{l) ~ l~ l appears. In a scattering 
process (13.181) . g(l) contains a term proportional to Z~ 3 . In order to see the non-analytic behavior 
of these integrals, we divide the integral Q(M) into a regular part Q(M;L) and a (possibly) 
singular part Q(M; oo) as 

pL pea 

Q(M) = Q(M; L) + Q(M; oo) = / dl I(M,l)g(l) +/ dl I(M,l)g(l), (A.12) 

Jo Jl 

where L is chosen to be large enough so that one can make the approximation g(l) = al~ n . By 
changing the integration variable to s = ML, the second part becomes 

Q(M,oc) = M ds^-g{s/M)^aM n+1 ds (A.13) 

J ML s J ML S" - 



19 



The approximation in the second equality is valid even in the massless limit because s/M is 
always larger than L and g(s/M) ~ a(s/M)~ n is a good approximation. The M dependence 
only comes from the integration around s = ML. For a small M, one can expand J(s) ps 
s/2 — (s) 3 /16 + ■ • • near s = ML, and Q(M, oo) becomes 

_ -Q k ^2k— n 



Q(M, oo) « aM s+1 / rfs V 

^ ML 



2 2k+1 k\(k + l) 



k=0 

M 2/c +2 l _^ c j ML . (A.14) 

^ 2 2fc+1 fc!(/c + 1)! V2A; - n + 1 J 

It is a regular function for an even integer n. For an odd integer n, the logarithmic term log M 
arises. The coefficient of the logarithmic term log M is given by 

2 n(^) !( n±l), ^ ) 

It is determined by the behavior of g(l) ~ od~ n at infinity and the behavior of the Bessel 
function J\{s) near s = 0. As an example, due to the asymptotic behavior ( 13.101) . the integral 
( 13. 6p at a = 1 is reduced to the above integral Q{M) with n = 1 and a = — e 2 /27r, and the 
logarithmic dependence in (13.111) is reproduced. 

When g(l) is an oscillating function, the behavior of the integral becomes different. It is 
the case corresponding to the circular motion in (13.231) and (13.241) . In this case, g(l) contains 
oscillating factors like sin kl dbt I — > 00. Suppose that g(l) behaves like g(l) ~ I n sin kl for 
I — > 00. If we divide the integral Q(M) into Q(M,L) and Q(M, 00) as above, a non-analytic 
behavior, if exists, must come from the second part Q(M, 00). By changing the integration 
variable to s — ML, it becomes 



00 



Q(M, 00) - aM^ I ds Jl(g)5i " ( f g/M) . (A.16) 

J ML sH 

Since ks/M > kL ^> 1, the term sin (ks/M) is oscillating rapidly. Since Ji(s) ~ s for small s, 
M-dependence of Q(M, 00) is determined by the integration, 

f°° sin ks/M , „ ,-. f°° sin jfc* 

/ — rfs = M n+1 / — ^— dt, (A.17) 

which is a power function of M. Thus no non- analytic behavior appears when g(l) is an 
oscillating function at infinity. 

B Existence of a positive region of F(a, M) 

In this appendix we prove that, in the uniform acceleration, whatever /(/) we choose the 
backreaction must take a positive value for a sufficiently large g = a 2 /AM 2 . This can be 
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shown as follows. If we take a different choice of f(l) from the simplest choice f(l) = 1/2, the 
backreaction changes as in (13.131) . This change of the backreaction must satisfy the condition 



lim SF(M) = (B.18) 

since it should vanish in the large M limit. Hence, for any small e, there exists M e that satisfies 
\5F(M)\ < e for M > M e . On the other hand, as shown in Figure. (H the function F(g) is 
positive for g > 1 and grows as a function of g and there exists g e such that F{g) > e for g > g e . 
For a fixed M > M e , we can always take a sufficiently large so that g = a 2 /AM 2 > g e . Then in 
such values of a and M, the backreaction becomes positive; 



F(a h M L 



<5F(M,) + F(-y >0. (B.19) 



e 2 a l M l v l/ MM; 
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